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ABSTRACT: We investigate implications of decoherence for quantum 
systems which are classically chaotic. We show that, in open systems, the 
rate of von Neumann entropy production quickly reaches an asymptotic 
value which is: (i) independent of the system-environment coupling, (ii) 
dictated by the dynamics of the system, and (hi) dominated by the largest 
Lyapunov exponent. These results shed a new light on the correspondence 
between quantum and classical dynamics as well as on the origins of the 
"arrow of time." 
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The relation between classical and quantum chaos has been always somewhat 
unclear 1 and, at times, even strained 2 . The cause of the difficulties can be traced 
to the fact that the defining characteristic of classical chaos - sensitive dependence 
on initial conditions - has no quantum counterpart: it is defined through the be- 
havior of neighboring trajectories 3 , a concept which is essentially alien to quantum 
mechanics. Moreover, when the natural language of quantum mechanics of closed 
systems is adopted, an analogue of the exponential divergence cannot be found. This 
is not to deny that many interesting insights into quantum mechanics have been 
arrived at by studying quantized versions of classically chaotic systems 4 . These 
insights have typically much to do with the energy spectra, and leave the issue of 
the relationship between the quantum and the classical largely open. 

The aim of this paper is to investigate implications of the process of decoherence 
for quantum chaos. Decoherence is caused by the loss of phase coherence between 
the set of preferred quantum states in the Hilbert space of the system due to the 
interaction with the environment 5 . Preferred states are singled out by their stabil- 
ity (measured, for example, by the rate of predictability loss - the rate of entropy 
increase) under the joint influence of the environment and the self-hamiltonian *\ 
Thus, the strength and nature of the coupling with the environment play a crucial 
role in selecting preferred states, which - given the distance-dependent nature of 
typical interactions - explains the special function of the position observable 
Coupling with the environment also sets the decoherence timescale - the time on 
which quantum interference between preferred states disappear [5- 9] . Classicality 
is then an emergent property of an open quantum system. It is caused by the inces- 
sant monitoring by the environment, the state of which keeps a "running record" 
of the preferred observables of the evolving quantum system. For simple quantum 
systems the programme sketched above can be carried out rigorously, and yields 
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intuitively appealing results 5 ' . For example, preferred states of an underdamped 
harmonic oscillator turn out to be its coherent states 8 . 

If decoherence does induce a transition from quantum to classical, then it should 
be possible to utilize it in the context of quantum chaos to establish a more straight- 
forward correspondence between the behavior of classically chaotic systems and 
their quantum counterparts. With this goal in mind we will consider a classically 
chaotic system, characterized by a potential V(x), coupled to an external environ- 
ment. A master equation for the density operator of an open quantum system can 
be derived under a variety of reasonable assumptions 10 . Here, we shall focus our 
considerations on the simplest special case, the high temperature limit of an ohmic 
environment, leaving the discussion of other cases (low temperature regime, non- 
ohmic environments, etc) for later publications. In this case, the Wigner function 
of the system evolves according to 10 : 

W = {H,W} PB + J2 22 n J + ; d* n + 1 Vd% n+1 W + 2 1 d p (pW) +Dd 2 pp W (1) 
n>l 

where 7 is the relaxation rate, the diffusion coefficient is D = 2^mk B T (T is the 
temperature of the environment and m is the mass of the system) . The first term is 
the Poisson bracket, which generates the ordinary Liouville flow. Both the Poisson 
bracket and the higher derivative terms result from an expansion of the Moyal 
bracket, {H,W}mb = ~ * sm(ih{H, W}p B )/h, which generates evolution in phase 
space of a closed system (this expansion is valid when V(x) is an analytic function). 
The last two terms in (1) arise due to the interaction with the environment. The 
first of them produces relaxation - gradual loss of energy to the reservoir - and the 
last one diffusion (this diffusion term is responsible for the decoherence process) . 

As is clear from (1), as a consequence of the quantum correction terms that 
contain higher derivatives, the Wigner function of an isolated non-linear system 
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does not follow the classical Liouville flow. In a classically chaotic system these 
non-classical corrections rapidly gain importance, as can be seen by the follow- 
ing argument: When a chaotic flow is investigated locally in the phase space, the 
evolution operator can be expanded in coordinates "comoving" with a reference 
trajectory. The pattern of flow of the neighboring trajectories is then generated by 
the Jacobian of the transformation. Eigenvalues of this Jacobian are known as local 
Lyapunov exponents Aj, which must sum to zero, since the transformation preserves 
phase space volume. Eigenvectors define directions in the phase space along which 
the neighboring trajectories either only expand (A^ > 0) or only contract (A^ < 0) 
with respect to the fiducial trajectory at a rate given by the corresponding Lyapunov 
exponent 3 . The exponential contraction rapidly generates small scale structure in 
the Wigner function. Thus, the high derivative terms are d p W oc a p n W with 
a p oc <J P (0) exp(At) where A is a Lyapunov exponent. Hence, non classical correc- 
tions will become important after a characteristic crossover time t x which can be 
estimated by comparing the magnitude of the nonlinear corrections with the contri- 
bution of the Poisson bracket in equation (1). Defining a characteristic length for 

the non-linear terms in the potential as \ n oc (d x V / d^: +1 V) l l n we obtain that the 

(n) 

n-th order term in (1) becomes comparable with the Poisson bracket at a time 
given by 

4 n) ocA- 1 log( X n^(0)/n) (2) 

One of the points we want to make here is that, in the presence of decoherence, 
the regime in which the quantum correction become important is easily avoided. 
The main reason for this is the existence of diffusive effects which put a lower bound 
on the small scale structure which can be produced by the chaotic evolution. As 
a result, o~ p (t) can never become sufficiently small to result in large corrections: 
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Poisson bracket is an excellent approximation of the quantum Moyal bracket for 
smooth Wigner functions. This line of reasoning, as we will argue below, has also 
important consequences concerning the rate of entropy production. 

To simplify our analysis, based on equation (1), we will neglect the relaxation 
term which, as pointed out in the literature 9 , can be made arbitrarily small without 
decreasing the effectiveness of the decoherence process (e.g. by letting 7 approach 
zero while keeping D constant). In this way, we will focus in the important reversible 
classical limit 6 ' 9 > n . We will not limit ourselves to models leading to Eq. (1) 
which, as they break the symmetry between x and p coupling with the environment 
through position, have momentum diffusion only. It is convenient (especially in 
the context of quantum optics, where the "rotating wave approximation" can be 
invoked) to use a symmetric coupling (of the form a)b + ab^ , where a and b are the 
annihilation operators of the system and the mode of the environment 12 ). The 
corresponding equation differs from (1) in the form of the diffusion which is now 
symmetric, oc D(dp p + d% x )W. We shall alternate between using this symmetric 
diffusion and the more exact diffusion operator of equation (1) in the discussion 
below. 

Our objective here is to study the interplay between the evolution which clas- 
sically results in an exponential divergence of neighboring trajectories (the char- 
acteristic feature of chaos) and the destruction of quantum coherence between a 
preferred set of states in the Hilbert space of an open system (a defining feature 
of decoherence). We shall do this by using a simple unstable system that still cap- 
tures the essential features we want to consider. In general, a chaotic geodesic flow 
pattern is locally analogous to the one occurring near a saddle point, with stable 
and unstable directions defined in an obvious manner. The simplest example of 
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such a saddle point is afforded by an unstable harmonic oscillator. We shall use it 
as a "generic" model of locally chaotic phase space dynamics in our considerations 
below. In this case, the potential is V(x) = (A is the Lyapunov exponent) 

and equation (1) reduces to a very simple form. We will use the lessons learned from 
this simple unstable system to argue that quantum corrections, which in the case 
of the unstable oscillator vanish identically, can be neglected whenever decoherence 
is effective. Our analysis will also show that three different stages of the evolution 
can be identified: (i) decoherence, (ii) approximately reversible Liouville flow and 
(iii) irreversible diffusion-dominated evolution. 

To analyze in detail the unstable oscillator it is convenient to use contracting 
and expanding coordinates defined by " = (p =F rnXx). Evolution generated by Eq. 
(1) causes exponential expansion in v and, without the diffusive term, it would also 
cause an exponential contraction in u, so that the volume in the phase space (as 
well as entropy) would be constant (see Fig. 1). Expansion in v would also result 
in an exponential decrease of gradients in that direction. Thus, after a sufficient 
number of e-foldings the equation governing evolution of W would be dominated 
by the expression: 



The characteristic dispersion, which will play an important role below, is 



We can now easily examine the fate of a generic initial quantum state. The 
general solution of equation (3) can be found by noticing that the eigenfunctions 
of the operator appearing in its right hand side are v n F m (u/a c ) where F m (x) = 
exp(— /2)H m -\(x / \/2) and H m (x) are Hermite polynomials. Expanding W(u, v, t) 




(3) 
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in terms of these eigenfunctions (whose eigenvalue is simply —(n + m)) we obtain 

W(u, v,t) = J2 a nm (ve- Xt ) n F m (u)e~ mujt (5) 

m>l 

From this expression we see that the Wigner function depends on v only through 
the combination vq = vexp(— At), which is the comoving coordinate. That is, along 
this direction, the Wigner function just expands. Moreover, after a few dynamical 
times the most important contribution to (5) will always come from the m = 1 term. 
Thus, in the contracting direction the Wigner function approaches a Gaussian with 
a critical width and has the form: 

1 2 In 2 r°° 

W{u, v, t) « e~ u /2(Tc e~ w< / duW(u, v , t = 0). (6) 

y2na^ J—oo 

The existence of the critical width a c is a consequence of the interplay between 
the exponential divergence of trajectories and diffusion. In effect, a competition 
between the chaotic evolution (which attempts to "squeeze" the wavepacket in the 
contracting direction) , and the diffusion (which has the opposite tendency) leads to 
a compromise steady state which results in the Gaussian written above. 

How do these general considerations imply the existence of the three distin- 
guishable phases of the evolution? The analysis of decoherence, responsible for the 
quantum to classical transition, follows simply. Non-classical states possessing a 
rapidly oscillating non-positive W quickly evolve towards a mixture of localized 
states eventually resulting in a positive Wigner function. For example, if the initial 
state is a superposition of two coherent states separated by a distance L (along 
u), the ratio between the wavelength of the interference fringes i oc 1/L and a c is 
<Jc/£ 2 = (^/X)(L/X ( i^) 2 where is the thermal de Broglie wavelength 9 . There- 
fore, the decoherence time is 

r*e = 7- 1 (^f) 2 (7) 
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which, for macroscopic scales, is much smaller than the dynamical times even for 
very weakly dissipative systems. Our analysis of the decoherence period is still 
incomplete since equation (6) suggests that the negativity of the Wigner function 
may persist along the expanding direction (in that equation the initial state is not 
changed along v but just "stretched" by the geodesic flow). However, equation (3) 
was obtained by neglecting gradients along v. When a diffusion term Ddy V W is 
added to the right hand side of (3), the eigenfunctions change in such a way that 
the powers (ve~^ t ) n are replaced by Hermite polynomials H n (v/s/2a c )e~ n ^ t . If we 
reexpress the solution in terms of v q we notice that, for times of the order of 1/ A only 
the highest power in H n (v/y/2a c ) survives. This implies that the asymptotic form 
of W is no longer given by (6), which just contains the initial state expanded along 
the unstable direction. The correct expression simply contains a smoothed version 
of the initial condition in which the details smaller than a c are washed out along 
the unstable direction. Again, oscillations with wavelength £ oc 1/L are typically 
destroyed after the decoherence time (7). 

The analysis of the reversible and irreversible stages can be illustrated by fol- 
lowing the evolution of a Gaussian W. Here, the existence of a c is again very 
important. The von Neumann entropy H of a Gaussian state can be easily related 
to the area A enclosed by a 1-cr contour of the Wigner function. Thus, H is a 
monotonic function of A which, when A ^> h = 2nh approaches Tt ~ k^hiA/h. 
Using the above equations, one can show that the rate of entropy production is 

where cr p (t) is the width of the Gaussian along the direction of p. Consequently 
when the width of the Gaussian approaches the critical value, the entropy growth 
becomes equal to the positive Lyapunov exponent A. The evolution of H can be 
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approximately analyzed as follows: one can use (1) to show that the ratio R = 

(cjp(t)/o- c ) 2 evolves according to the equation R = 2A(1 — (3R) where (3(t), which is 

related to the degree of squeezing of the state, approaches unity exponentially fast. 

Solving this equation approximately (using (3 = 1) one gets: 

/ (7^(0) x-1 
H = A(l + (^2--l)exp(-2At)J (9) 
<j c 

Evolution of a classical distribution (corresponding to a non-negative W, ini- 
tially smooth on a scale much larger than a c and spread over a regular patch with 
A^> h) will typically proceed in two different stages (see Fig. 2). The first stage 
will be approximately area-preserving with the evolution dominated by the Liouville 
operator. It will last as long as each of the dimensions of the patch is much larger 
than the critical width. During this stage diffusion does little to alter the form of 
W. The Wigner function is merely "stretched" or "contracted" by the geodesic flow 
so that, with respect to the co-moving coordinates, "nothing happens" to W. By 
contrast, when the dimension of the patch becomes comparable with a c , diffusion 
will begin to dominate. Further contraction will be halted at a c but the stretching 
will proceed at the rate set by the positive Lyapunov exponent. As a result, the 
area (or, more generally, the volume) in phase space will increase at the rate set by 
(9) with <jp = a c . Using our approximate equation (9) one can estimate the time 
corresponding to the transition from reversible to irreversible evolution: 

r c = A- 1 ln(^l) (10) 

An important condition must be valid in order to apply the above arguments, 
illustrated here in the simplest example, to more complicated non-linear systems: 
we need to assure that the Wigner function follows approximately the evolution 
generated by the Poisson bracket (i.e. that the effect of the higher derivative terms 
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in (1) is small). As we pointed out above, in the absence of decoherence the quantum 
corrections to the evolution of W become important at a crossover time t x given 
by equation (2). These quantum corrections will remain small if diffusion is strong 
enough to prevent the formation of small structure in W. As we argued above, the 
formation of small structure induced by the exponential contraction is stopped after 
a time r c . Therefore, the condition for the Wigner function to evolve classically is 



which, taking into account equations (2) and (10), can be rewritten in the following 
suggestive way: 



This condition - a key criterion to assure the correspondence between quantum and 
classical dynamics - assures that W follows the Liouville flow (albeit with diffusive 
contributions) and allows us to apply the conclusions of our previous analysis. 

We have demonstrated that chaotic quantum systems can exhibit, in addition to 
the very rapid onset of decoherence, a nearly reversible phase of evolution which is 
necessarily followed by an irreversible stage in which the entropy increases linearly 
at the rate determined by the Lyapunov exponents. By contrast, open quantum sys- 
tems with regular classical analogs continue to evolve with little entropy production 
(although possibly with a significant change in dynamics 13 ). This nearly constant 
rate of (von Neumann) entropy production, a consequence of the interplay between 
the chaotic dynamics of the system and its interaction with the environment, sug- 
gests not only a clear distinction between the integrable and chaotic systems, but 
also shows that increase of entropy in the context of quantum measurement 14 and 
the dynamical aspects of the second law are intimately related and can be traced 



t x > t c , 



(11) 



\ n O c > h. 



(12) 
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to the same cause: Impossibility of isolating macroscopic systems from their envi- 
ronments. 

REFERENCES 

1. B.V. Chirikov, Phys. Rep. 52, 263 (1979); G. M. Zaslavsky, Phys. Rep. 80, 
157 (1981); M-J. Giannoni, A. Voros, J. Zinn-Justin, eds., Chaos and Quantum 
Physics, Les Houches Lectures LII (North Holland, Amsterdam, 1991). 

2. J. Ford and G. Mantica, Am. J. Phys. 60, 1086 (1992). 

3. See, e.g. G. loos, R. H. G. Hellmann and R. Stora eds., Chaotic Behavior in 
Deterministic Systems, Les Houches Lectures XXXVI (North Holland, Amster- 
dam, 1983); R. S. Mc Kay and J. O. Meiss eds., Hamiltonian Dynamical Systems 
(Hilger, Philadelphia, 1987). 

4. M. B. Berry, Proc. Roy. Soc. London, A423, 219 (1989); M. Gutzwiller, Chaos in 
Classical and Quantum Mechanics (Springer Verlag, New York, 1990); F. Haake, 
Quantum Signature of Chaos (Springer Verlag, New York, 1990) and references 
therein. 

5. W. H. Zurek, Phys. Rev. D24, 1516 (1981); ibid, D26, 1862 (1982); E. Joos and 
H. D. Zeh, Zeits. Phys. B59, 229 (1985). 

6. W. H. Zurek, Physics Today 44, 36 (1991); ibid, 46, 81 (1993); Prog. Theor. 
Phys. 89, 281 (1993). 

7. J. P. Paz, S. Habib and W. H. Zurek, Phys. Rev. D47, 488 (1993). 

8. W. H. Zurek, S. Habib and J. P. Paz, Phys. Rev. Lett. 70, 1187 (1993). 

9. W. H. Zurek, in Frontiers of Nonequilibrium Statistical Mechanics, G. T. Moore 
and M. O. Scully eds., (Plenum, New York, 1986). 

11 



10. A. O. Caldeira and A. J. Leggett, Physica 121 A, 587 (1983); W. G. Unruh and 
W. H. Zurek, Phys Rev. D40, 1071 (1989); B. L. Hu, J. P. Paz and Y. Zhang, 
Phys. Rev. D45, 2843 (1992); ibid~D47, 1576 (1993). 

11. M. Gell-Mann and J. B. Hartle, Phys. Rev. D47, 3345 (1993). 

12. W. H. Louisell, Quantum Statistical Properties of Radiation (Wiley, New York, 
1973); V. Buzek and P. Knight; Quantum interference, superposition states of 
light and nonclassical effects, Prog, in Opt. (1993), to appear. 

13. G.J. Milburn and C.A. Holmes, Phys. Rev. Lett. 56, 2237 (1986). 

14. J. von Neumann, Mathematical Foundations of Quantum Mechanics, english 
translation by R. T. Beyer (Princeton Univ. Press, Princeton, 1955); H. D. Zeh, 
The Direction of Time, (Springer- Verlag, New York, 1991). 



12 



Figure Captions 

Figure 1: Three stages of the evolution of an initially Gaussian Wigner function 
shown in two different coordinate systems. Fig. (a)-(c) display W in physical 
coordinates (u,v, with v in logarithmic scale). For an open system which evolves 
according to equation (1) with V = Ax 2 /2 (forefront), the width of the distribution 
reaches the asymptotic value a c . By contrast, when W evolves unitarily (densely 
hatched and shown in the back), it continues to be squeezed in u. Fig. (a')-(c') 
shows the same three stages of the evolution, but now in co-moving coordinates 
(v = vexp(\t),u = uexp(-Xt)). In the unitary case (cross-hatched, shown in 
the back) W does not change. The interaction with the environment causes an 
exponential increase in the apparent width of W (forefront). Since the width of the 
Gaussians in the expanding direction is approximately the same, the asymptotic 
regime of the diffusive evolution leads to an exponential increase of the area enclosed 
in l-o" contour. Consequently, the entropy increases linearly at a rate determined 
by the Lyapunov exponent. 

Figure 2: The rate of von Neumann entropy production for the quantum open 
system. The initial state is a Gaussian for which H(t = 0) 3> 1 and with the initial 
width along the contracting direction much larger than a c . The (nearly) reversible 
and irreversible stages of the evolution are clearly distinguished. 
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